In silico screening of Pueraria tuberosa (PTY-2) for targeting COVID-19 by countering dual targets Mpro and TMPRSS2.

COVID-19 pandemic was started in Wuhan city of China in December 2019; immensely affected global population. Herein, an effort was made to identify potential inhibitors from active phytochemicals of Pueraria tuberosa (PTY-2) via molecular docking study. Our study showed five potential inhibitors (Robinin, Genistin, Daidzin, Hydroxytuberosone, Tuberostan) against Mpro and five inhibitors (Robinin, Anhydrotuberosin, Daidzin, Hydroxytuberosone, Stigmasterol) against TMPRSS2. Out of these, Robinin, Daidzin and Hydroxytuberosone were common inhibitors for Mpro and TMPRSS2. Among these, Robinin showed the highest binding affinity, therefore, tested for MD simulation runs and found stable. ADMET analysis revealed the best-docked compounds are safe and follow the Lipinski Rule of Five. Thus, it could be suggested that phytochemicals of PTY-2 could serve as potential inhibitors for COVID-19 targets.Communicated by Ramaswamy H. Sarma.

One microsecond MD simulations of the SARS-CoV-2 main protease and hydroxychloroquine complex reveal the intricate nature of binding.

Currently, several vaccines and antivirals across the globe are in clinical trials. Hydroxychloroquine (HCQ) was reported to inhibit the SARS-CoV-2 virus in antiviral assays. Here, it raises the curiosity about the molecular target of HCQ inside the cell. It may inhibit some of the viral targets, or some other complex mechanisms must be at disposal towards action mechanisms. In some of the viruses, proteases are experimentally reported to be a potential target of HCQ. However, no in-depth investigations are available in the literature yet. Henceforth, we have carried out extensive, one-microsecond long molecular dynamics simulations of the bound complex of hydroxychloroquine with main protease (Mpro) of SARS-CoV-2. Our analysis found that HCQ binds within the catalytic pocket of Mpro and remains stable upto one-third of simulation time but further causes increased fluctuations in simulation parameters. In the end, the HCQ does not possess any pre-formed hydrogen bond, other non-covalent interactions with Mpro, ultimately showing the unsteadiness in binding at catalytic binding pocket and may suggest that HCQ may not inhibit the Mpro. In the future, this study would require experimental validation on enzyme assays against Mpro, and that may be the final say. Communicated by Ramaswamy H. Sarma.

Amyloidogenic proteins in the SARS-CoV and SARS-CoV-2 proteomes.

The phenomenon of protein aggregation is associated with a wide range of human diseases. Our knowledge of the aggregation behaviour of viral proteins, however, is still rather limited. Here, we investigated this behaviour in the SARS-CoV and SARS-CoV-2 proteomes. An initial analysis using a panel of sequence-based predictors suggested the presence of multiple aggregation-prone regions (APRs) in these proteomes and revealed a strong aggregation propensity in some SARS-CoV-2 proteins. We then studied the in vitro aggregation of predicted aggregation-prone SARS-CoV and SARS-CoV-2 proteins and protein regions, including the signal sequence peptide and fusion peptides 1 and 2 of the spike protein, a peptide from the NSP6 protein, and the ORF10 and NSP11 proteins. Our results show that these peptides and proteins can form amyloid aggregates. We used circular dichroism spectroscopy to reveal the presence of ß-sheet rich cores in aggregates and X-ray diffraction and Raman spectroscopy to confirm the formation of amyloid structures. Furthermore, we demonstrated that SARS-CoV-2 NSP11 aggregates are toxic to mammalian cell cultures. These results motivate further studies about the possible role of aggregation of SARS proteins in protein misfolding diseases and other human conditions.

An insight into SARS-CoV-2 membrane protein interaction with spike, envelope, and nucleocapsid proteins.

Intraviral protein-protein interactions are crucial for replication, pathogenicity, and viral assembly. Among these, virus assembly is a critical step as it regulates the arrangements of viral structural proteins and helps in the encapsulation of genomic material. SARS-CoV-2 structural proteins play an essential role in the self-rearrangement, RNA encapsulation, and mature virus particle formation. In SARS-CoV, the membrane protein interacts with the envelope and spike protein in Endoplasmic Reticulum Golgi Intermediate Complex (ERGIC) to form an assembly in the lipid bilayer, followed by membrane-ribonucleoprotein (nucleocapsid) interaction. In this study, we tried to understand the interaction of membrane protein's interaction with envelope, spike, and nucleocapsid proteins using protein-protein docking. Further, simulation studies were performed up to 100 ns to examine the stability of protein-protein complexes of Membrane-Envelope, Membrane-Spike, and Membrane-Nucleocapsid proteins. Prime MM-GBSA showed high binding energy calculations for the simulated structures than the docked complex. The interactions identified in our study will be of great importance, as it provides valuable insight into the protein-protein complex, which could be the potential drug targets for future studies.Communicated by Ramaswamy H. Sarma.

Mitoxantrone dihydrochloride, an FDA approved drug, binds with SARS-CoV-2 NSP1 C-terminal† † Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ra07434b

One of the major virulence factors of SARS-CoV-2, NSP1, is a vital drug target due to its role in host immune evasion through multiple pathways. NSP1 protein is associated with inhibiting host mRNA translation by binding to the small subunit of ribosome through its C-terminal region. Previously, we have shown the structural dynamics of the NSP1 C-terminal region (NSP1-CTR) in different physiological environments. So, it would be very interesting to investigate the druggable compounds that could bind with NSP1-CTR. Here, in this article, we have performed different spectroscopic technique-based binding assays of an anticancer drug mitoxantrone dihydrochloride (MTX) against the NSP1-CTR. We have also performed molecular dynamics simulations of the docked complex with two different force fields up to one microsecond. Overall, our results have suggested good binding between NSP1-CTR and MTX and may have implications in developing therapeutic strategies targeting the NSP1 protein of SARS-CoV-2. SARS-CoV-2 NSP1 is believed to be a potential drug target. Its C-terminal region which associates with ribosomal subunit has shown good binding with an anticancer drug, Mitoxantrone.

Mitoxantrone dihydrochloride, an FDA approved drug, binds with SARS-CoV-2 NSP1 C-terminal.

One of the major virulence factors of SARS-CoV-2, NSP1, is a vital drug target due to its role in host immune evasion through multiple pathways. NSP1 protein is associated with inhibiting host mRNA translation by binding to the small subunit of ribosome through its C-terminal region. Previously, we have shown the structural dynamics of the NSP1 C-terminal region (NSP1-CTR) in different physiological environments. So, it would be very interesting to investigate the druggable compounds that could bind with NSP1-CTR. Here, in this article, we have performed different spectroscopic technique-based binding assays of an anticancer drug mitoxantrone dihydrochloride (MTX) against the NSP1-CTR. We have also performed molecular dynamics simulations of the docked complex with two different force fields up to one microsecond. Overall, our results have suggested good binding between NSP1-CTR and MTX and may have implications in developing therapeutic strategies targeting the NSP1 protein of SARS-CoV-2.

Microsecond simulations and CD spectroscopy reveals the intrinsically disordered nature of SARS-CoV-2 spike-C-terminal cytoplasmic tail (residues 1242-1273) in isolation.

All available SARS-CoV-2 spike protein crystal and cryo-EM structures have shown missing electron densities for cytosolic C-terminal regions (CTR). Generally, the missing electron densities point towards the intrinsically disordered nature of the protein region (IDPR). This curiosity has led us to investigate the cytosolic CTR of the spike glycoprotein of SARS-CoV-2 in isolation. The spike CTR is supposed to be from 1235 to 1273 residues or 1242-1273 residues based on our used prediction. Therefore, we have demonstrated the structural conformation of cytosolic region and its dynamics through computer simulations up to microsecond timescale using OPLS and CHARMM forcefields. The simulations have revealed the unstructured conformation of cytosolic region. Further, we have validated our computational observations with circular dichroism (CD) spectroscopy-based experiments and found its signature spectra at 198 nm. We believe that our findings will surely help in understanding the structure-function relationship of the spike protein's cytosolic region.

Investigating the conformational dynamics of SARS-CoV-2 NSP6 protein with emphasis on non-transmembrane 91-112 & 231-290 regions.

The NSP6 protein of SARS-CoV-2 is a transmembrane protein, with some regions lying outside the membrane. Besides a brief role of NSP6 in autophagosome formation, this is not studied significantly. Also, there is no structural information available to date. Based on the prediction by TMHMM server for transmembrane prediction, it is found that the N-terminal residues (1-11), middle region residues (91-112), and C-terminal residues (231-290) lies outside the membrane. Molecular Dynamics (MD) simulations showed that NSP6 consists of helical structures. In contrast, the membrane outside lying region (91-112) showed partial helicity, which was further used as a model and obtained disordered type conformation during 1.5 µs. Additionally, a 200ns simulation study of residues 231-290 have shown significant conformational changes. As compared to helical and beta-sheet conformations in its structure model, the 200ns simulation resulted in the loss of beta-sheet structures while helical regions remained intact. Further, we have experimentally characterized the residue 91-112 by using reductionist approaches. CD spectroscopy suggests that the NSP6 (91-112) is disordered-like region in isolation, which gains helical conformation in different biological mimic environmental conditions. These studies can be helpful to study NSP6 (91-112) interactions with host proteins, where different protein conformations might play a significant role. The present study adds up more information about the NSP6 protein aspect, which could be exploited for its host protein interaction and pathogenesis.

SARS-CoV-2 NSP1 C-terminal (residues 131-180) is an intrinsically disordered region in isolation.

The NSP1- C terminal structure in complex with ribosome using cryo-EM is available now, and the N-terminal region structure in isolation is also deciphered in literature. However, as a reductionist approach, the conformation of NSP1- C terminal region (NSP1-CTR; amino acids 131-180) has not been studied in isolation. We found that NSP1-CTR conformation is disordered in an aqueous solution. Further, we examined the conformational propensity towards alpha-helical structure using trifluoroethanol, we observed induction of helical structure conformation using CD spectroscopy. Additionally, in SDS, NSP1-CTR shows a conformational change from disordered to ordered, possibly gaining alpha-helix in part. But in the presence of neutral lipid DOPC, a slight change in conformation is observed, which implies the possible role of hydrophobic interaction and electrostatic interaction on the conformational changes of NSP1. Fluorescence-based studies have shown a blue shift and fluorescence quenching in the presence of SDS, TFE, and lipid vesicles. In agreement with these results, fluorescence lifetime and fluorescence anisotropy decay suggest a change in conformational dynamics. The zeta potential studies further validated that the conformational dynamics are primarily because of hydrophobic interaction. These experimental studies were complemented through Molecular Dynamics (MD) simulations, which have shown a good correlation and testifies our experiments. We believe that the intrinsically disordered nature of the NSP1-CTR will have implications for enhanced molecular recognition feature properties of this IDR, which may add disorder to order transition and disorder-based binding promiscuity with its interacting proteins.

Naturally Occurring Bioactives as Antivirals: Emphasis on Coronavirus Infection.

The current coronavirus disease (COVID-19) outbreak is a significant threat to human health and the worldwide economy. Coronaviruses cause a variety of diseases, such as pneumonia-like upper respiratory tract illnesses, gastroenteritis, encephalitis, multiple organ failure involving lungs and kidneys which might cause death. Since the pandemic started there have been more than 107 million COVID-19 infections caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and ∼2.4 million deaths globally. SARS-CoV-2 is easily transmitted from person-to-person and has spread quickly across all continents. With the continued increase in morbidity and mortality caused by COVID-19, and the damage to the global economy, there is an urgent need for effective prevention and treatment strategies. The advent of safe and effective vaccines has been a significant step forward in the battle against COVID-19, however treatment of the symptoms associated with the disease still requires new anti-viral and anti-inflammatory drug therapies. To this end, scientists have been investigating available natural products that may be effective against SARS-CoV-2, with some products showing promise in fighting several viral infections. Since many natural products are dietary components or are prepared as dietary supplements people tend to consider them safer than synthetic drugs. For example, Traditional Chinese Medicines have been effectively utilized to treat SARS-CoV-2 infected patients with promising results. In this review, we summarize the current knowledge of COVID-19 therapies and the therapeutic potential of medicinal plant extracts and natural compounds for the treatment of several viral infections, with special emphasis on SARS-CoV-2 infection. Realistic strategies that can be employed for the effective use of bioactive compounds for anti-SARS-CoV-2 research are also provided.

Conformational dynamics of 13 amino acids long NSP11 of SARS-CoV-2 under membrane mimetics and different solvent conditions.

The intrinsically disordered proteins/regions (IDPs/IDPRs) are known to be responsible for multiple cellular processes and are associated with many chronic diseases. In viruses, the existence of a disordered proteome is also proven and is related to its conformational dynamics inside the host. The SARS-CoV-2 has a large proteome, in which, structure and functions of all proteins are not known yet, along with non-structural protein 11 (nsp11). In this study, we have performed extensive experimentation on nsp11. Our results based on the CD spectroscopy gives characteristic disordered spectrum for IDPs. Further, we investigated the conformational behavior of nsp11 in the presence of membrane mimetic environment, α-helix inducer, and natural osmolyte. In the presence of negatively charged and neutral liposomes, nsp11 remains disordered. However, with SDS micelle, it adopted an α-helical conformation, suggesting the helical propensity of nsp11. Finally, we again confirmed the IDP behavior of nsp11 using MD simulations. In future, this conformational dynamic study could help to clarify its functional importance in SARS-CoV-2 infection.

Reprofiling of approved drugs against SARS-CoV-2 main protease: an in-silico study.

Given the COVID-19 pandemic, currently, there are many drugs in clinical trials against this virus. Among the excellent drug targets of SARS-CoV-2 are its proteases (Nsp3 and Nsp5) that plays vital role in polyprotein processing giving rise to functional nonstructural proteins, essential for viral replication and survival. Nsp5 (also known as Mpro) hydrolyzes replicase polyprotein (1ab) at eleven different sites. For targeting Mpro, we have employed drug repurposing approach to identify potential inhibitors of SARS-CoV-2 in a shorter time span. Screening of approved drugs through docking reveals Hyaluronic acid and Acarbose among the top hits which are showing strong interactions with catalytic site residues of Mpro. We have also performed docking of drugs Lopinavir, Ribavirin, and Azithromycin on SARS-CoV-2 Mpro. Further, binding of these compounds (Hyaluronic acid, Acarbose, and Lopinavir) is validated by extensive molecular dynamics simulation of 500 ns where these drugs show stable binding with Mpro. We believe that the high-affinity binding of these compounds will help in designing novel strategies for structure-based drug discovery against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Translation-Associated Mutational U-Pressure in the First ORF of SARS-CoV-2 and Other Coronaviruses.

Within 4 months of the ongoing COVID-19 pandemic caused by SARS-CoV-2, more than 250 nucleotide mutations have been detected in ORF1ab of the virus isolated from infected persons from different parts of the globe. These observations open up an obvious question about the rate and direction of mutational pressure for further vaccine and therapeutics designing. In this study, we did a comparative analysis of ORF1a and ORF1b by using the first isolate (Wuhan strain) as the parent sequence. We observed that most of the nucleotide mutations are C to U transitions. The rate of synonymous C to U transitions is significantly higher than the rate of non-synonymous ones, indicating negative selection on amino acid substitutions. Further, trends in nucleotide usage bias have been investigated in 49 coronaviruses species. A strong bias in nucleotide usage in fourfold degenerate sites toward uracil residues is seen in ORF1ab of all the studied coronaviruses: both in the ORF1a and in the ORF1b translated thanks to the programmed ribosomal frameshifting that has an efficiency of 14 - 45% in different species. A more substantial mutational U-pressure is observed in ORF1a than in ORF1b perhaps because ORF1a is translated more frequently than ORF1b. Mutational U-pressure is there even in ORFs that are not translated from genomic RNA plus strands, but the bias is weaker than in ORF1ab. Unlike other nucleotide mutations, mutational U-pressure caused by cytosine deamination, mostly occurring during the RNA plus strand replication and also translation, cannot be corrected by the proof-reading machinery of coronaviruses. The knowledge generated on the mutational U-pressure that becomes stronger during translation of viral RNA plus strands has implications for vaccine and nucleoside analog development for treating COVID-19 and other coronavirus infections.

Exploring the SARS-CoV-2 structural proteins for multi-epitope vaccine development: an in-silico approach.

INTRODUCTION: The ongoing life-threatening pandemic of coronavirus disease 2019 (COVID-19) has extensively affected the world. During this global health crisis, it is fundamentally crucial to find strategies to combat SARS-CoV-2. Despite several efforts in this direction and continuing clinical trials, no vaccine has been approved for it yet. METHODS: To find a preventive measure, we have computationally designed a multi-epitopic subunit vaccine using immuno-informatic approaches. RESULTS: The structural proteins of SARS-CoV-2 involved in its survival and pathogenicity were used to predict antigenic epitopes. The antigenic epitopes were capable of eliciting a strong humoral as well as cell-mediated immune response, our predictions suggest. The final vaccine was constructed by joining the all epitopes with specific linkers and to enhance their stability and immunogenicity. The physicochemical property of the vaccine was assessed. The vaccine 3D structure prediction and validation were done and docked with the human TLR-3 receptor. Furthermore, molecular dynamics simulations of the vaccine-TLR-3 receptor complex are employed to assess its dynamic motions and binding stability in-silico. CONCLUSION: Based on this study, we strongly suggest synthesizing this vaccine, which further can be tested in-vitro and in-vivo to check its potency in a cure for COVID-19.

Understanding COVID-19 via comparative analysis of dark proteomes of SARS-CoV-2, human SARS and bat SARS-like coronaviruses.

The recently emerged coronavirus designated as SARS-CoV-2 (also known as 2019 novel coronavirus (2019-nCoV) or Wuhan coronavirus) is a causative agent of coronavirus disease 2019 (COVID-19), which is rapidly spreading throughout the world now. More than 1.21 million cases of SARS-CoV-2 infection and more than 67,000 COVID-19-associated mortalities have been reported worldwide till the writing of this article, and these numbers are increasing every passing hour. The World Health Organization (WHO) has declared the SARS-CoV-2 spread as a global public health emergency and admitted COVID-19 as a pandemic now. Multiple sequence alignment data correlated with the already published reports on SARS-CoV-2 evolution indicated that this virus is closely related to the bat severe acute respiratory syndrome-like coronavirus (bat SARS-like CoV) and the well-studied human SARS coronavirus (SARS-CoV). The disordered regions in viral proteins are associated with the viral infectivity and pathogenicity. Therefore, in this study, we have exploited a set of complementary computational approaches to examine the dark proteomes of SARS-CoV-2, bat SARS-like, and human SARS CoVs by analysing the prevalence of intrinsic disorder in their proteins. According to our findings, SARS-CoV-2 proteome contains very significant levels of structural order. In fact, except for nucleocapsid, Nsp8, and ORF6, the vast majority of SARS-CoV-2 proteins are mostly ordered proteins containing less intrinsically disordered protein regions (IDPRs). However, IDPRs found in SARS-CoV-2 proteins are functionally important. For example, cleavage sites in its replicase 1ab polyprotein are found to be highly disordered, and almost all SARS-CoV-2 proteins contains molecular recognition features (MoRFs), which are intrinsic disorder-based protein-protein interaction sites that are commonly utilized by proteins for interaction with specific partners. The results of our extensive investigation of the dark side of SARS-CoV-2 proteome will have important implications in understanding the structural and non-structural biology of SARS or SARS-like coronaviruses.